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Abstract. We report on a device capable of imaging second-order spatio-temporal 
correlations g^'^'^ (x, r) between photons. The imager is based on a monolithic array of single- 
photon avalanche diodes (SPADs) implemented in CMOS technology and a simple algorithm 
to treat multiphoton time-of-arrival distributions from different SPAD pairs. It is capable of 80 
ps temporal resolution with fluxes as low as 10 photons/s at room temperature. An important 
application might be the local imaging of g^'^'^ as a means of confirming the presence of true 
Bose-Einstein macroscopic coherence (BEC) of cavity exciton polaritons. 
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Figure 1. Interference fringes (a) at A = 770 nm wavelength used to verify the BEC of 
polaritons in GaAs microcavity (reproduced from Ref. |4|) and (b) at A = 546 nm measured 
for the green Hne of pulsed Hg-Ar discharge lamp. 



Recent experiments HI El El ID have reported the Bose-Einstein condensation (BEC) 
phase transition in polariton system in a semiconductor microcavity. The macroscopic 
quantum degeneracy is typically detected by probing the statistical properties of light emitted 
from a microcavity, under the presumption that the statistics of the exciton polaritons are 
faithfully transferred to the emanating photons. 

It has been further assumed that observation of the interference fringes similar to those 
in Michelson or Young interferometers (figure [Ita)) is sufficient to establish the fact of 
macroscopic coherence in polariton system. [[B SI E]] Two points on the wave front separated 
by a distance Xi2 produce intensity pattern I1+I2 + 2\/ Iil2\g^^^ {^12 , ^) | cos(A99i2), such that 
the fringe visibility measures the magnitude of the first-order correlation function. But simply 
measuring this quantity alone is ambiguous because a coherent light source {e.g., a photon 
laser or decaying polariton BEC) can exhibit the same first-order correlations as a chaotic (or 
thermal) light source {e.g. Hg-Ar discharge lamp in figure [B^b)). Table [H shows that proper 
disambiguation of a coherent state also requires measurement of the second-order correlation 
function g^'^\i^i2,T)=^^'^^^^^^^ associated with intensity noise correlations. Here, /i,2(i) 
is the light intensity at a point ±|xi2 and time t. The minimal condition to confirm the 
BEC phase transition in a polariton system then becomes 5'^^^(x, 0) = 5'^^^(x, 0) = 1 (third 
column of Table [I]). So far, very few theoretical and experimental studies of the second-order 
correlations in a polariton system were limited to the k=Q point in the momentum space (in 
the lateral cavity direction), reporting thus a spatially-averaged value of ^'^^^(0) and ignoring 
the fact that 2D BEC can be achieved only in spatially confined systems. As a consequence, 
the model of Ref.[6J predicts increasing correlation peak height g'^^^ (0) with polariton number 
caused by the strong scattering effects above the critical threshold density. Such behaviour 
for polaritons at /c = has been confirmed in Ref. [7 J, while the experimental results of 
Ref. [8 J disagree with such behaviour. To confirm the BEC phase transition in a polariton 
system, one needs to distinguish the presence of both coherent condensate (^(^)=^(^)=1) and 
thermal noncondensate (s'^^Ms'^^^) fractions, which can be achieved by local measurement of 
the first-order and second-order correlations. However, the small dimensions of the polariton 
distribution in the microcavity, means that such measurements must be capable of resolving 
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Table 1. Values of first and second order correlation functions for incoherent, coherent and 
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a spatial dependence in ^'(^^ All these requirements demand an integrated monolithic photon 
detector, like a camera, but capable of imaging intensity noise correlations. 

In this letter, we present an imaging device allowing the second-order (intensity) 
correlations of a field to be probed locally. The current implementation comprises a 4 x 4 
silicon array of single-photon avalanche diodes (SPADs) implemented in 0.35 /xm CMOS 
technology [9J. The device incorporates on-chip high bandwidth I/O circuitry (figure 
[2]) for off-chip data processing of multi photon arrivals. We demonstrate the operational 
performance of our monolithic ^^^^^ -imager in a miniature replication of the Hanbury Brown 
and Twiss (HBT) stellar interferometer [lOJ. 

Each pixel in our SPAD array is based on an avalanche photodiode operating in the 
Geiger mode. The 3.5/xm-diameter active region of a SPAD n-well pixel consists of a 
p+-n junction reverse-biased above its breakdown voltage. When a photon is absorbed 
in the multiplication region, an avalanche is triggered with a certain probability. The 
measured single photon detection probability of such 0.35/xm-CMOS SPAD pixels is 25 
% at 546 nm wavelength and 4 V excess voltage above the breakdown threshold. (The 
detection probability is 40% at the peak sensitivity wavelength 450 nm [9J.) The avalanche 
breakdown is subsequently quenched by an on-chip ballast resistor, which is used to read 
out the photodetection events. Its value defines the dead time of the detector which is 
15 ns for each SPAD pixel of the 4x4 array considered here. The chip comprises built- 
in high-bandwidth electronics to convert Geiger pulses into digital signals for off-chip data 
processing. This design drastically improves the signal-to-noise ratio. At room temperature, 
the lowest detectable photon flux is set by the dark count rate (DCR) of SPADs in the 5- 
10 Hz range. Such low DCR is achieved by using small n-wells of diameter 3.5 /xm. The 
lowest detectable photon flux density in our experiments is ^10^ photons / s-cm'^ , well below 
the level in polariton BEG experiments. The array pitches are 30 jum horizontal and 43 fim 

vertical. All 16 detectors in the array have separate parallel outputs so that ^ ^ = 120 

simultaneous pairwise measurements are possible at a temporal resolution limited by the 
SPAD jitter characteristics (80 ps). 

The pairwise intensity noise correlations g^'^'^ (x^j , r) are computed using an external four- 
channel 6 GHz bandwidth digital oscilloscope (Wavemaster 8600A, LeCroy) by programming 
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Figure 2. Exploded micrograph of the 4 x 4 SPAD array, 
it with an algorithmic version of equation ([T]): 

M f 

NMY, x['^\n)Xf'\n+l) 

9^'\^3^r)= ? (1) 

E E ^^^(^) E Xf'\n'+l) 

where integers enumerate detector pixels, Xi and Xj are discrete random variables 

whose values (no event) or 1 (photon detection) correspond to the binary data stream (figure 
WL^)) emanating from any pair of detectors Di and Dj, respectively. The spatial lag xij is set 
by the separation of the detector pair within the SPAD array. Time-lag increments r = IT are 
set by multiples of temporal resolution T, where NT is the width of measurement window 
and M is the overall number of measurements series. 

Unlike conventional detection methods based on start-stop timing histograms of 
delayed single photon arrivals [figure [S^a)], our approach implements properly normalized 
multiphoton distribution which is robust against missing detection events, the impact of 
Poisson-like distribution decay ocexp(— /xr) at large r and intensity modulation. Equation 
([T]) permits any count rates and temporal window of interest and does not require a statistical 
hypothesis to normalize g^'^\ Figure [3tb) shows the benefits of our technique by comparing 
g^'^^r) measurements of incoherent broad-band light using multiphoton arrivals and the 
standard two-photon histogram. 

In figure [S^b), for the standard approach based on time-delayed coincidence events, the 
histogram was acquired and then normalized assuming that at the lag r=0 there should 
be no correlations for incoherent light source (i.e., assuming that g^'^\0)=l). We shall 
remind that a probability distribution of time intervals between two consecutive photons is 
p{0\t) = /xexp(— /xr) for a Poisson process [[HI [111 and that for a multi-mode Gaussian 
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Figure 3. Comparison of single-photon and multi-photon measurements, (a) Traces from two 
detectors at 5 MHz count rate. Standard histogram of delayed single-photon arrivals shows a 
large number of lost detections, (b) Computed g^'^'^ based on multi-photon arrivals (green) and 
the standard technique (red) showing the impact of lost detections at count rates /i ~2.5 MHz 
at each detector and the width of measurement window NT =100ns. 



(chaotic light) and a coherent (Poisson statistics) states these distributions are the same [13]. 
Therefore, when normahzed in such way, the standard method shows the onset of Poisson 
distribution g^'^^r) = p(0|r)//i 1—^1^1 for large time lag |r| ^ Thus in fig 

Ob), the average count rate at detectors is /x = 2.5 MHz, and for r = 50 ns, the error is 
A^'^^) ^ 0. 13. The conventional procedure is thus limited by the small count rate and temporal 
window width product. At the same time, our approach produces the correct correlations 
independent of the measurement interval NT and the photon flux intensity. 

The s'^^) -imager was tested by measuring the statistical properties of an extended thermal 
light source. As a model system for a quasi-monochromatic chaotic light source, we used a 
Hg-Ar spectral wavelength calibration lamp (bulb CAL-2000-B, Ocean Optics) with U-folded 
discharge and cold cathode. This lamp, designed for operation in the AC regime, also operated 
well with a DC power supply (160 V @ 15 mA). To start the discharge, we used the original 
AC power supply of CAL2000 source which was connected in parallel with a DC source via a 
filter (2H inductance) such that the AC supply was gradually turned off while the DC source 
was gradually turned on. 

The light emitted by the lamp was transmitted through a 10 nm bandpass filter (FL543.5- 
10, Thorlabs), which keeps only the emission at the green line of mercury (546 nm). The light 
then was injected into aim long multimode fiber (figureffia)) with core diameter w=200 iim. 
The other end of the fiber was used to illuminate the SPAD array placed in the far field zone of 
the fiber end, at a distance L=2 cm from the fiber [figure |4l^a)]. The numerical aperture of the 
fiber (NA=0.22) assumes that the whole 4x4 SPAD array is over illuminated. Such extended 



Quantum Imager 



(3) ^Hg-Ar 



SPAD array 




multimode fiber 
(200Mm core) 



nm 



filter 



CJ 



00 



DO 
D5 
D9 
D15 



160V (DC) 



2cm 



(b) 2 



-g 





(C) 2. 



(2) 



'''''''' 

1 - A^ - ^V-A^u^^-VV^ 



,(2) 



=' D0&D15 
1 - />^ ^ /-^^^^^r^^y^.,7^,N^/V^/^ 



Time, ns 



— 1 — 

-40 



-20 



-r- 






20 

D1 



D5 



40 

D9 



D13 



50 100 150 200 
Detector separation x, [nm] 



D13 



Jm 



Figure 4. (a) SPAD array configured as miniature HBT interferometer for measuring 200 fim 
(b) correlations between adjacent detectors (x^^q = 30/im) and (c) extreme diagonal detectors 
(^0,15 = 158/im). (d) Measured (points) and calculated (curve) second-order in function of 
detector separation The micrograph inset shows the position of detectors within the array, (e) 
imaged second-order correlation maxima along the row of array in (d). It is assumed that 
^(2) =2 along the diagonal. 



thermal light source is of the angular width w / L=10"^ rad and exhibits first-order correlations 
g^^"^ (figure [H^b)), when SPAD array is replaced by Young double-pinhole interferometer (not 
shown in the figure). 
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The second-order spatio-temporal correlations for a non-polarized single-mode chaotic 
light source are determined by the first-order correlations [[HI [Ml [15]] with coherence time 
rc=2\/27r In2/A6j due to the inhomogeneous broadening Acj (FWHM) of the emission line 

where the second term in the right hand side takes into account the decorrelation effects 
due to unpolarized light (the coefficient 1/2), the zero-delay degree of spatial coherence and 
the Gaussian profile of the delayed first-order correlation function for an inhomogeneously- 
broadened line. 

To examine the temporal correlations, we used detector pairs in the middle of the array 
as well pairs at diagonally opposite corners, thus providing correlations g^'^^Xij.r) between 
two different regions separated by 30 and 158 /xm, respectively. At small separation (figure 
12b)), the data shows an excess of coincidences g^'^'^ (^5^9, 0)=1.25. At large distance, detector 
counts are uncorrected (figure 12c)). 

At temporal resolution of the scope T=l ns, the measured coherence time Tc is 5.2 ns 
such that the impact of integration effects on measured correlation peak width and height is 
less than 2 %. fW\ The corresponding linewidth of the source is 130 MHz (FWHM), which 
is less than the Doppler width of the green line of Hg and can be attributed to the Dicke 
linewidth narrowing due to the buffering effect of Ar in the lamp bulb. 

The spatial oscillations of g^'^\xij,0) due to the spatial coherence factor in ([2]) were 
measured by selecting detector pairs from a row of the array (D1-D13 row in figure [2]) with 
detector separation 30, 60 and 90 /xm. Figure |4][d) shows the correlation excess at zero-delay 
is well fitted by a sine function in equation (O, yielding the angular width of the source 
0.9x 10~^ rad, close to the estimated value. 

Being limited by the number of acquisition channels, we were able to record 
simultaneous correlations between four independent detectors. In figure |4](e), 5'(^^(x^j,0) 
measured along the array row is plotted as a pairwise correlation map g^'^'^ (i, j). In this image 
map, the spatial oscillations of the coherence factor are clearly visible. 

In conclusion, we have presented a ^^^^^ -imager built with conventional CMOS 
technology, which is capable of measuring second-order spatio-temporal correlated photons 
and thereby offers an important means for verifying the existence of a BEC state of cavity 
exciton polaritons. Future work will include the development of larger arrays of SPADs, the 
integration of on-chip data processing based on equation ([T]), and the extension to other g^'^'>- 
imaging applications. 
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